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ABSTRACT 

We present Chandra data from a 3 1 .7 ks observation of the Centaurus cluster, using the ACIS- 
S detector. Images of the X-ray emission show a plume-like feature at the centre of the clus- 
ter, of extent 60 arcsec (20 kpc in projection). The feature has the same metallicity as gas at 
a similar radius, but is cooler. Using adaptive binning, we generate temperature, abundance 
and absorption maps of the cluster core. The radial abundance profile shows that the previ- 
ously known, steep abundance gradient peaks with a metallicity of 1 .3 — 1 .8 Z p at a radius of 
about 45 arcsec (15 kpc), before falling back to 0.4 Z s at the centre of the cluster. A radial 
temperature profile shows that the temperature decreases inwards. We determine the spatial 
distributions of each of two temperature components, where applicable. The radiative cooling 
time of the cooler component within the inner 10 arcsec (3 kpc) is less than 2 x 10 7 yr. X-ray 
holes in the image coincident with the radio lobes are seen, as well as two outer sharp tem- 
perature drops, or cold fronts. The origin of the plume is unclear. The existence of the strong 
abundance gradient is a strong constraint on extensive convection or gas motion driven by a 
central radio source. 

Key words: X-rays: galaxies — galaxies: clusters: individual: Centaurus — intergalactic 
medium — cooling flows 



1 INTRODUCTION 

The Centaurus cluster, Abell 3526, at a redshift z = 0.0104 (Lucey, 
Currie & Dickens 1986) is a nearby, X-ray bright galaxy cluster. 
Analysis of ASCA, and earlier data (e.g. Allen et al. 2001) indi- 
cates it hosts a cooling flow (see Fabian 1994 for a review) of 
30-50 M yr _1 , centred on the cD galaxy NGC 4696. Previous 
observations of Centaurus (e.g. Allen & Fabian 1994) show a sys- 
tem with smooth, elliptical, X-ray isophotes, indicating the system 
is relatively relaxed. However, there is evidence for a current or 
past merger event (Allen & Fabian 1994; Churazov et al. 1999; Fu- 
rusho et al. 2001). The neighbouring subcluster, Cen 45, centred 
on NGC 4709, has a velocity which is 1500 km s -1 higher than the 
main Centaurus cluster, also known as Cen 30 (Lucey, Currie & 
Dickens 1986). 

Observations of the cluster using ROSAT and ASCA show that 
the central region of the cluster is rich in metals, with a large abun- 
dance gradient (Fukazawa et al. 1994; Ikebe et al. 1998; Allen et 
al. 2001). 
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2 THE X-RAY IMAGE 

The data presented here are from a 31.7 ks observation of the Cen- 
taurus cluster by the Chandra Observatory, using the ACIS-S de- 
tector. The temperature of the focal plane was — 120°C at the time 
of the observation. No period in the observation was affected by 
flares, so no data were excluded from the analysis. The interaction 
of the radio source PKS 1246-410 with the X-ray plasma is dis- 
cussed elsewhere (Taylor, Fabian & Allen 2001). All errors quoted 
are l-o" uncertainties. Distances assume H Q = 50 km s _1 Mpc -1 . 
Data were fitted using the XSPEC package (version 1 1 .0. 1 ; Arnaud 
1996) and metallicities are relative to solar abundances defined by 
Anders & Grevesse (1989). 

Fig. |l| shows a coloured raw (~ 1 arcsec pixels) image of the 
inner 140 arcsec of the S3 chip. Events are coloured according to 
the energy recorded (0.5-1 keV red, 1-2 keV green, 2-7 keV blue). 
The plume-like structure to the NE resembles a twisted sheet, and 
was seen in an observation of the cluster using the ROSAT HRI 
(Sparks, Jedrzejewski & Macchetto 1993). 

Fig. |^ shows an adaptively-smoothed (Ebeling, White & Ran- 
garajan 2001) image of the raw data. It has been smoothed with a 
minimum significance of 3-<J. The image highlights that the emis- 
sion is diminished from the north-east to the south-west close to 
the nucleus. It clearly shows the spiral plume-like feature, appear- 
ing to emerge from the nucleus in the south, and heading towards 
the north-east. The plume is approximately 60 arcsec long (20 kpc 
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Figure 1. Colour X-ray image of the inner part of the ACIS-S3 chip centred 
on the core of the Centaurus cluster. Events between 0.5 and 1 keV are 
coloured red, those between 1 and 2 keV green, and those between 2 and 
7 keV blue. The pixel size is ~ 1 arcsec. The image measures ~ 140 arcsec 
(46 kpc) vertically. The emission peaks at around 110 counts per square- 
arcsecond between 0.5 and 7 keV. 



in projection). The energy-coloured image (Fig. |lj) shows the centre 
of the cluster, and the plume, is cooler (redder) than the rest of the 
cluster. The two deficits or 'holes' in emission immediately to the 
NE and NW of the core are associated with the lobes of the radio 
source (Taylor et al. 2001). The continuation of the holes to the SE 
and SW is probably also related with the southern extension to the 
radio lobes. 

Fabian et al. (2001b) hypothesise that the filament-like feature 
seen in X-ray emission in the core of the Abell 1795 cluster may 
be caused by the relative motion of the cD galaxy to the rest of 
the cluster (estimated to be at least 150 kms ). The results of 
Dickens, Currie & Lucey (1986) show that NGC 4696 lies close to 
the mean velocity of Cen 30. However, the radio maps of Taylor et 
al. (2001) show that the radio lobes bend to the south of the core. 
This indicates that NGC 4696 is moving north relative to the local 
intracluster medium. The aspect of the radio lobes could include a 
projection effect, with a velocity component along our line of sight 
of 100-200 kms- 1 . 

Observations of the core of the cluster at optical wavelengths 
show a dust-lane in NGC 4696 (Shobbrook 1963). It has been 
noted that the X-ray emission has a similar morphology to the dust 
lane (Sparks et al. 1993). Using four HST archived observations 
(datasets U3560101T, U3560102T, U62G8401R and U62G8402R) 
of the galaxy at 5407 A and 7940 A, we have made a ratio image of 
the optical emission in two bands highlighting the dust-lane. This 
is shown in Fig. ^, overlayed with a contour map of the centre of 
the adaptively smoothed image, Fig. ^ The swirling morphology 
of the dust-lane is similar to, but smaller than, the X-ray plume. 

We find no evidence for hard X-ray emission (between 2-7 or 
4-7 keV) from the core of the radio source PKS 1246-410 (posi- 
tion taken from Taylor et al. 2001) at the 2-a level. Hard X-ray 
emission was seen in Centaurus using ASCA (Allen, Di Matteo & 
Fabian 2000), which has a hundredfold larger point-spread function 
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Figure 2. Adaptively smoothed image of the centre of the Centaurus cluster 
in the energy range 0.5 to 7 keV with ~ 0.5 arcsec pixels. The image mea- 
sures 175 arcsec vertically, which corresponds to a distance of 58 kpc. The 
minimum ASMOOTH significance is 3-0". 
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Figure 3. Optical HST image showing the dust lane in the centre of 
NGC 4696 overlayed by contours from the adaptively smoothed image of 
the X-ray data (Fig. The optical image shows the ratio of the emission 
in a band with a central wavelength of 5407 A relative to one at 7940 A. 

than Chandra. There is tentative evidence for diffuse hard X-ray 
emission in the Chandra data which will be studied elsewhere. 

3 SPECTRAL FITS 

3.1 Spectrum of central 10 arcsec 

We extracted a spectrum from the central 10 arcsec radius of 
the cluster. The spectrum was fitted with an absorbed single- 
temperature (1-T) MEKAL model (Mewe, Gronenschild & van den 
Oord 1985; Liedahl, Osterheld & Goldstein 1995), in the energy 
range of 0.5 to 7 keV, where the Chandra response-matrix is well 
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calibrated. The spectrum was binned to ensure a minimum number 
of counts per spectral bin, chosen here to be 22 arbitrarily. We gen- 
erated an effective response and ancillary matrix, averaged from the 
matrices applicable to each part of the region, but weighted accord- 
ing to the relative number of counts. The background was removed 
using a blank field background-spectrum of the same part of the 
detector, created using the prescription of Markevitch (2001). Al- 
though the background fields have lower galactic absorption (/V H = 
0.02 x 10 22 cm" 2 ) than the region observed (0.081 x 10 22 cm -2 , 
obtained using the FTOOLS NH program; Dickey & Lockman 1990) 
its level is insignificant compared with the bright cluster emission 
studied here. 

The results of the fit are summarised in Table |l|. We find a 
temperature of 0.71 keV, a metallicity of O.16Z , absorbed by 
N H — 0.26 x 10 22 cm~ 2 . The 1-T model poorly fits the spectrum 
(Reduced x 2 = 660/128 = 5.2). 

In regions of cooling flow clusters where radiative cooling 
dominates, i.e. where the radiative cooling time is shorter than 
the age of the flow (about 3-5 Gyr; Allen & Fabian 1997, corre- 
sponding to about one arcmin radius in the Centaurus cluster), the 
X-ray surface brightness profiles indicated, if there is no balanc- 
ing heat source, that the inferred mass deposition of cooled gas is 
roughly proportional to radius (M oc r). This requires that the gas 
there is multiphase (Nulsen 1986; Thomas, Fabian & Nulsen 1987). 
The spectrum of multiphase gas should be fit better with a two- 
temperature (2-T) model than a 1-T model (e.g. Allen et al. 2001). 
An even better model should be a cooling flow (CF: Johnstone et 
al. 1992), possibly with intrinsic absorption. 

Studies of several cooling flow clusters using the XMM- 
Newton RGS (e.g. Peterson et al. 2001; Tamura et al. 2001) have 
shown that soft X-ray emission lines expected from gas cooling to 
low temperatures are absent in the data. A model (cooling flow to 
a minimum temperature, CFMT) which does fit the data, is a cool- 
ing flow in which gas cools from the upper temperature to a lower 
temperature, which may be one third to one fifth of the upper one. 
Such a model requires that the gas disappears once it reaches the 
minimum temperature, perhaps because of intrinsic absorption for 
which the absorbed CF model is an approximation, because of mix- 
ing or because of rapid heating which turns the cool phase into the 
hot one (see Fabian et al. 2001a). 

In this paper we try all these models (1-T, 2-T, CF and CFMT) 
on the data. Previous work with ASCA has shown that more than 
one temperature component is required in the core of the Centau- 
rus cluster (Fukazawa et al. 1994; Ikebe et al. 1998; Allen et al. 
2001). Projection effects due to hotter outer gas being along the 
line of sight to cooler gas are also examined for the plume. We do 
not however carry out an exhaustive set of modelling due to the 
unknown 3-dimensional geometry of the cooler gas and the limita- 
tions of the models themselves, which do not include effects due for 
example to gravitational heating in a flow. Moreover, as we show, 
the metallicity of the gas is spatially complex, which further com- 
plicates detailed modelling. It is our intention here to identify the 
gross behaviour of the gas in terms of its temperature and abun- 
dance distributions with position. 

We fitted the spectrum of the core with a 2-T MEKAL model, 
fixing the abundances of the two phases to be the same as there is 
not enough information in the spectrum to measure the abundances 
separately and accurately. The results of this fit are also shown in 
Table [l]. The temperatures of the two phases are around 0.7 keV 
and 1.4 keV. The metallicity (~ 0.6 Z e ) is much larger than for the 
1-T result (~ 0.2 Z Q ). The quality of the fit is considerably better 
than the 1-T modei (Reduced % 2 = 200/132 = 2.4), but still not 



good. Allowing the two abundances to vary did not significantly 
improve the fit. The best-fitting metallicities were Z (cold- 

phase) and 0.46+gjjg Z,.~ (hot-phase), but the temperatures were un- 
affected. Adding a third temperature to the 2-T model did not sig- 
nificantly improve the fit. However, adding a third component fixed 
at 4 keV to account for projected emission was an improvement. 
The radiative cooling time of the cooler component in the 2-T fit is 
about 2 x 10 7 yr. 

The fits can be improved further by switching from a MEKAL 
model to an APEC model (Brickhouse et al. 2000). The reduced x 2 
for an APEC 2-T fit is 2.0. The residuals of the fit show there is some 
excess emission above 3 keV. This emission can be fit either using 
a power-law or a high-temperature plasma. We show in Table |l| the 
slight changes induced by adding a plasma of fixed temperature 
4 keV and metallicity 0.4 Z . Adding this component is justified 
as we expect contributing emission from the cluster in front of the 
core. 

The spectrum can also be fitted (but with a poorer goodness- 
of-fit) using an intrinsically absorbed cooling flow model (RJCOOL; 
Johnstone et al. 1992), plus a plasma fixed at 4 keV and 0.4 Z Q . 
We absorbed both components by a fixed galactic absorber, and al- 
lowed the intrinsic absorption of the cooling flow to vary. Normally 
we would add another plasma component to account for the local 
emission of the cluster near the cooling gas, and gravitational work 
done. However, at the centre of the cluster the cooling flow domi- 
nates, so no extra component was necessary. 

We derive a cooling flow mass deposition rate of 1 1 M p yr~' . 
The gas cools from a temperature of 1.1 keV and has a metal- 
licity of 0.7 Z Q . The cooling flow is intrinsically absorbed by 
2.4 x 10 21 cm . The spectrum, fitted with the cooling flow model, 
is shown in Fig. Q The cooling flow model calculates its emission 
using MEKAL plasmas. We speculate that an APEC version of the 
code would fit the spectra better, but unfortunately this code is still 
under development. 

Using a CFMT model (Table |l|) we find the mass deposition 
rate and upper temperature to be unchanged, a decrease in the in- 
trinsic absorption of the cooling flow to a value similar to that found 
using 2-T fits, and a lower cooling temperature of 0.40 keV. The fit 
is better than a single cooling flow, but not as good as a 2-T model. 

All the fits we have tried to the core spectrum show a large 
residual around 0.75 keV, independent of the model used (shown 
for the RJCOOL model in the lower panel of Fig. [J). MEKAL and 
APEC models with varying elemental abundances do not improve 
the fit to this feature. The emission is well fitted with a 5-function. 
By examining the spectra of annuli, we have determined that this 
feature is not significant outside a radius of 10 arcsec (3 kpc) from 
the centre. It is also present in the inner 5 arcsec. The effective 
area of the telescope varies smoothly over the energy range of the 
emission, making the feature unlikely to be a calibration error. 

There are also negative-residuals just above 2 keV. The effec- 
tive area per unit energy of the telescope almost halves in the range 
2 to 2.1 keV, so this feature may be a calibration issue. Accounting 
for these two spectral features would greatly improve the fit to the 
data, and would lead to a much smaller reduced % 2 . 

We show in Fig. ^ the contributions towards x 2 for three dif- 
ferent fits to the spectrum of the core. The feature at 0.75 keV, in 
all spectra, has been removed by fitting it with a Gaussian. We also 
included a plasma fixed a 4 keV and 0.4 Z r , l to model the projected 
cluster emission. The top panel shows a 2-T MEKAL model, the 
middle panel a 2-T APEC model, and the lower panel an intrinsi- 
cally absorbed RJCOOL model with fixed galactic absorption. 
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Figure 4. Spectrum of the inner 1 arcseconds of the cluster, fitted with a 1 - 
T MEKAL model (fixed at 4 keV and 0.4 Z „), plus an intrinsically absorbed 
RJCOOL cooling flow. 



3.2 Spectrum of the plume 

We extracted a spectrum from the plume-like region in the cluster. 
It, like the core, is poorly fit by a 1-T model. Using a 2-T model, 
the fitted temperature of the hot component is ~ 0. 1 keV larger than 
the hot component in the core. The cooler component has a similar 
temperature. The abundance of this region is significantly higher 
than the core (an increase of 0.2 — 0.3 Z Q ). 

We also fitted the spectrum of an area at a similar radius to the 
plume (but not including the plume) with a 2-T model. The results 
are shown in Table |^ as 'plume radius'. The plume appears to be 
~ 0.25 keV cooler (in both fitted temperatures). There is no signif- 
icant abundance difference. A cooling flow model (with a compo- 
nent to account for local cluster emission) gives a mass deposition 
rate of 4M ffi yr~' , with a high metallicity of 1.6Z Q . 

Our next task was to test whether the multi-temperature na- 
ture of the plume is a projection effect. We compared the emission 
measure of the upper temperature component of a 2-T model fit- 
ted to the plume, with the emission measure of a 1-T fit to the 
spectrum in another region with the same size and radius as the 
plume, but at an angle of 90° to it. The emission measure of the 
upper-temperature component of the plume was 70 percent greater 
than the emission measure of this comparison region. Therefore, 
the multi-temperature nature of the plume is not wholly a projec- 
tion effect but is intrinsic to the plume. 

By assuming a cylindrical geometry for the plume (20 kpc 
long, radius 2.3 kpc), and fitting a 1-T model to its spectrum, the 
emission measure implies an electron density of ~ 0. 13 cm~ 3 . The 
total gas mass of the plume is therefore ~ 10 9 M Q . This mass is 
compatible with the hypothesis that this gas was stripped off an in- 
falling galaxy. The mean cooling time of the gas in the plume is 
approximately 10 8 yr, so the gas, or galaxy, must have a velocity 
of > 200 km s in order to not have cooled. The similarity of the 
abundance of the plume to the other gas at that radius raises a prob- 
lem for the stripping hypothesis. 



4 IMAGING ANALYSIS 
4.1 Radial profiles 

In order to quantify the properties of the X-ray emitting gas as a 
function of radius, we analysed the spectra of circular annuli around 
the core of the cluster. This assumes that the properties of the gas 




Figure 5. Contribution to x 2 (times sign of the difference between model 
and spectrum) of fitting the spectrum of the inner 10 arcsec with three dif- 
ferent models, including the contribution of the emission at 0.75 keV. (top) 
2-T MEKAL model plus a fixed component at 4 keV and 0.4 Z~ . (middle) 2- 
T APEC model with a fixed component, (bottom) RJCOOL component plus 
another fixed component. 
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Figure 6. Spectrum of the plume, fitted with a 2-T APEC model. 



are largely symmetric, an assumption we will test in §4.2 and §[0|. 
Spectra were extracted, and appropriate background spectra, re- 
sponse and ancillary matrices were created as in § 3. 1 for every an- 



nulus. We fitted each spectrum with 1-T and 2-T MEKAL models, 
leaving metallicity and absorption as free parameters. To identify 
whether a 2-T fit was substantially better than a 1-T fit for a par- 
ticular annulus, we computed the F-statistic. If the probability of 
the observed value of F (or a larger value) was less than or equal to 
0.01, we chose the 2-T fit for that annulus. The annuli were allowed 
to vary in width to ensure that there were 8000 counts in each spec- 
trum (except for the outer annulus). Point sources were excluded 
from the analysis. 
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Area 



Counts Model 



(10 22 c 



kT, (keV) 



! ) 



kT 2 (keV) 

or M ( M Q yr~ 



Z(Z„) 



Reduced x 



Central 10 arcsec 17835 



phabs(mekal) 0.26±0.01 0.71 ±0.01 

phabs(apec) 0.23 ±0.01 0.76 ±0.01 

phabs(mekal + mekal) 0.19±0.01 0.67±0.01 1.41±0.05 

phabs(apec + apec) 0.16±0.01 0.68±0.01 1.45±0.03 

phabs(apec+apec+[apec])* 0.16 ±0.01 0.66 ±0.01 1.33 ±0.03 

phabs(zphabs(rjcool) + 0.24±0.02 1.14±0.02 M=11.3±0.7 
[mekal]) 1 ' 

phabs(zphabs(rjcool 0.16±0.02 1.15±0.03 M=11.0±0.7 

- rjcool) + [mekal])* 0.40 ± 0.02 



0.16±0.01 
0.14±0.01 
59+ 07 
0.62 ±0.06 



0.64 



-o.t 



-0.05 

' 1 n no 



0.49+ 



660/128 = 
671/128 = 
299/126 = 
256/126 = 
243/125 = 
371/127 = 



5.2 
5.2 
2.4 
2.0 
1.9 
2.9 



339/126 = 2.7 



Plume 



14206 phabs(mekal) 0.18 ±0.01 1.14 ±0.01 

phabs(apec) 0.17 ±0.01 1.22 ±0.01 

phabs(mekal + mekal) 0.17 ±0.01 0.79 ±0.02 1.58 ±0.04 

phabs(apec + apec) 0.16±0.01 0.82±0.02 1.62±0.03 

phabs(mekal + [mekal] + 1 .56 ± 0.06 

zphabs(rjcool)) 5 0.28 ±0.02 M = 4.1 ±0.3 

phabs(mekal + [mekal] + L94 -oi8 

zphabs(rjcool-rjcool)) 0.11±0.03 0.54±0.03 M==3.5±0.3 



0.28±0.01 
0.30 ±0.02 
8Q+0.10 

u -<"-0.08 

0.81 ±0.08 



1.58 



-0.35 
-0.23 



583/138 = 
451/138 = 
278/136 = 
219/136 = 
289/136 = 



4.2 
3.3 
2.0 
1.6 
2.1 



1.15 ±0.15 254/135 == 1.9 



Plume radius 



24694 phabs(apec + apec) 



0.16 ±0.01 



1.05 ±0.03 1.89 ±0.04 



0.97±0.07 281/192 == 1.5 



Table 1. Results of spectral fits. Model names are shown in XSPEC notation. *Fit with a 2-T model, with an additional component fixed at 4 keV and 0.4 Zq. 
1 1ntrinsically-absorbed cooling flow plus a plasma component fixed at 4 keV and 0.4 Z s , These are both absorbed by a fixed absorber, set to the galactic value 
of absorption (0.081 x 10 22 cm -2 ). The temperature shown is the upper temperature of the cooling flow. N H is the intrinsic cooling flow absorption. *Like 
(t), but subtracting another cooling flow component with the same normalisation and abundance, but different upper temperature as the first cooling flow. 
This simulates a cooling flow cooling from an upper temperature to a lower temperature (shown on the second line). §Like (t), but with an additional MEKAL 
component to account for the local cluster emission. Its temperature is set to the upper temperature of the cooling flow, and its abundance is set to the same as 
the cooling flow. 



Fig. [7] shows the radial profiles of temperature, abundance, ab- 
sorption and the reduced y} of the fits. The values of the reduced 
X 2 show the fits are relatively poor in the centre. Two temperature 
points are given at the radii where the F-statistic gave reason to 
do so, showing the temperatures of the two temperature compo- 
nents. The vertical error bars show the l-o" errors to the fit for each 
variable, whilst the horizontal error bars mark the maximum and 
minimum radii of each annulus. 

The temperature in the outer regions of the cluster is 
~ 3.4 keV. Towards the centre of the cluster, inside a radius of 
~ 50 arcsec (17 kpc), a 2-T model is required to fit the spectra. The 
two temperatures drop to 1.35 and 0.66 keV at the centre (consis- 
tent with §3.1). The abundance rises from a value of 

1.3 Z 



0.4 Z in the 
at a radius of ~ 45 arcsec 



outer regions, to a maximum of <• 
(15 kpc). It then drops in the centre to ~ 0.4Z Q . The intrinsic ab- 
sorption rises from a value of ~ 10 21 cm~ 2 in the outer regions, to 
~ 2 x 10 21 cm~ 2 at the centre. The lower-panel shows the temper- 
ature results of fitting with 2-T and 1-T APEC models, rather than 
MEKAL models. The abundance and absorptions are very similar to 
the MEKAL results, so we do not plot them here. 

Using models with variable abundances (VMEKAL and 
VAPEC), we find the relative abundances near the core are more 
similar to those found in Type la supernovae rather than in Type 
II supernovae (defined by ratios given in Dupke & Arnaud 2001), 
tieing together those elements which are produced by the same 
process in each type of supernova. An even better fit is found 
by allowing the abundances to be free: we find abundances (rel- 
ative to solar), for an annulus 35 arcsec away from the core us- 
ing VAPEC, of O=0.5, Mg=0.8, Ne = 2 x lO" 6 , Si=1.3, S=1.5, 
Ar=1.5, Ca=2.3, Fe=l.l, and Ni=1.8. VMEKAL gives different 
results: O=0.7, Mg=0.2, Ne=0, Si=1.8, S=2.3, Ar=2.4, Ca=3.4, 



Fe=2.1, and Ni=2.1. The level of Mg depends in detail on how the 
Fe L emission is modelled. 

The volume filling fraction, /, of the cooler component in the 
2-T fits is given by 



EM 2 
EM, 



(1) 



where the emission measures and temperatures of the cooler and 
hotter components are EM l , T, and EM 2 ,T , respectively. This as- 
sumes thermal pressure balance between the two components. The 
radial profile of / (Fig. ^) shows that the cooler gas is detected 
when it exceeds about one per cent of the volume; it exceeds 10 per 
cent within the inner 10 arcsec. 

We also fitted an absorbed cooling flow model to each an- 
nulus. We used an intrinsically absorbed cooling flow (RJCOOL) 
plus a MEKAL model, both absorbed by fixed galactic absorption. 
The temperature and metallicity of the plasma model were fixed 
to the upper temperature and abundance of the RJCOOL model. 
Fig. ^ (top) shows a plot of the cumulative mass deposition rate 
as a function of radius: the combined mass deposition rate of all 
the annuli inside a particular radius. The plot rises from a value 
of ~ 7 M e yr~' in the central 6 arcsec, to ~ 50 M yr _1 at large 
radii (300 arcsec; 100 kpc), in agreement with previous measure- 
ments for the mass deposition rate (Allen et al. 2001). Also shown 
in Fig. ^ are the intrinsic absorption of the cooling flow, its upper 
temperature, its metallicity, and the reduced x 2 of eac h fit. 

As with the spectrum of the core of the cluster, we fitted each 
annulus with a CFMT model (the minimum temperature for each 
annulus rises from 0.4 keV at the centre to 1 keV at a radius of 1 ar- 
cmin). This model is only a significant improvement to the standard 
CF model inside a radius of around 6.5 arcsec from the core. 
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Figure 7. Results of fits taken of spectra taken from radial annuli. Errors are 
1-c A radius of 1 arcsec corresponds to 2.04 detector coordinates, or a dis- 
tance of 0.33 kpc. Radii at which there are two points have been fitted with 
a 2-T spectral model. The second-lowest panel shows the results of fitting 
with APEC models, rather than MEKAL models. The lowest panel shows the 
volume filling factor, /, of the cool component (where appropriate). 



We note that the level of absorption found is always higher 
than the expected galactic value. This may, in part, be an ACIS- 
S calibration issue. When fitting cooling flow spectra it may also 
characterize the lack of low temperature emission lines, such as are 
absent from XMM- NewtonfRGS spectra of cooling flows. We do 
not have the spectral resolution or sensitivity to detect absorption 
edges which would demonstrate the presence of intrinsic absorp- 
tion. 
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Figure 8. Cumulative mass deposition profile, calculated by fitting the mass 
deposition rate in annuli, and summing from the centre outwards. Also 
shown are the intrinsic absorption of the cooling flow, its upper tempera- 
ture, the abundance of the gas, and the reduced %' of the fit. The radiative 



cooling time of the gas at 300 arcsec is 
< 10 9 yr within about 60 arcsec. 



10 yr, decreasing inward to 



4.2 East-west anisotropy 

The east side of the cluster contains the plume-like feature. To in- 
vestigate whether there is significant difference between the east 
and west sides of the cluster, we analysed annuli cut in a north- 
south direction. Fig. ^ shows the two temperature and abundance 
profiles for the cluster. The annuli were chosen to contain 20000 
counts between 0.5 and 7 keV, and then were split into east and 
west parts. 1-T a nd 2 -T MEKAL models were fitted using the same 
procedure as in i 4.1 and an F-test was used to decide whether the 



2-T model was a significant improvement over the 1-T model. 

The plot shows significant differences in abundance between 
the two sides at radii of ~ 40-50 arcsec. There are also sharp out- 
ward rises in temperature in the east and west, at radii of ~ 70 and 
160 arcsec, respectively. We will discuss these further in the next 
section. Between these radii, there are temperature differences be- 
tween the two halves of the cluster. These radii match the region 
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Figure 9. Radial abundance and temperature profiles of the east and west 
of the cluster. 



in MEKAL and APEC radial plots where two temperature compo- 
nents are often preferred (Fig. 0). The anisotropy shows that more 
detailed analysis is required to map out the structure of the cluster. 



4.3 Fits to adaptively binned areas 

We adaptively binned (Sanders & Fabian 2001) the central 3 arcmin 
of the cluster, using a fractional error threshold of 0.03. Adaptive 
binning uses a bin size which changes to ensure a minimum number 
of counts in each bin. A threshold of 0.03 gives 1 100 counts per bin 
(except in certain cases). Fig. |l(](a) shows the adaptively binned 
image of the cluster. Point-sources were excluded from the analysis 
and are marked as white areas. The core of the cluster is shown to 
the south-east of the centre of the image. 

We extracted spectra from each bin, and generated appropri 



ate background spectra, response and ancillary matrices, as in § 1 
We then fitted each bin with 1-T and 2-T MEKAL models, finally 
using an F-test to decide whether the 2-T fit was a significant im- 
provement over the 1-T fit. The results for the two types of fit were 
then combined to create images showing the spatial distribution for 
each fitted parameter. Fig. llOlb) shows either the fitted temperature 
from the 1-T fit, or the upper temperature of the 2-T fit. Fig. [ic|(c) 
shows the lower temperature of the 2-T fit, where appropriate. Figs. 
|Io|(d) and (e) show the fitted abundance and absorption levels from 
the relevant fit. Fig. [u](f) shows the appropriate reduced x 2 value. 

The value of a parameter in the images is indicated by a colour 
scale. Subtracted point sources are marked by the maximum value. 
We plot the l-a errors to the fit as a pair of diagonal boxes in the 
centre of each bin. The top-right box shows the upper bound, and 
the bottom-left shows the lower bound. Masked-out point sources 
are shown as areas coloured with the maximum value. 

2-T models are required to fit the spectra in the inner 30 arc- 
sec. The 2-T temperature results in the centre (~ 0.7 and 1.2 keV) 
are cooler than the 1-T results in the outer regions (~ 3 keV). The 
abundance of the gas in the outer regions has a similar value to that 
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Figure 11. Mean radiative cooling time of the gas, estimated from the X-ray 
surface brightness (see Fabian et al. 2000). The mean cooling time is less 
than 5 x 10 s yr inside a radius of 45 arcsec, below 10 8 yr inside 10 arcsec 
(3 kpc) radius, and falls to 3 x 10 7 yr within 5 arcsec of the centre. 



in the core (~ 0.5 Z Q ), but there is a rise in the metallicity about 



1 arcmin away from the core to ~ 1.3 Z , as in { 4. 1 



The images show the properties of the emitting gas are not 
spherically symmetric. The temperature sharply increases at be- 
tween 70 and 160 arcsec radius, depending on position angle 
(Fig. ^). Such features have been seen in other clusters, and are 
termed cold fronts (Markevitch et al. 2000). The most striking front 
is to the north-east of the core in the Centaurus cluster, where the 
emission suddenly drops. There is an extensive cooler region to the 
west of the core (~ 2.4 keV; Fig. |l0^b)). This area has a higher 
abundance than an equivalent region to the east (Fig. [u](d)). In ad- 
dition, there is more absorption (Fig. |l(](e)) to the west. 

In the innermost core of the cluster (at a radius of less than 
4 arcsec) there are a few bins for which a 1-T fit is preferred using 
an F-test (probabilities for F of ~ 0.03). Owing to this, they have 
much lower abundances than the surrounding bins. If their spectra 
are combined, a 2-T fit is easily preferred, and the results of the fit 
closely match the values of the surrounding bins. This is likely to 
be due to the difficulties of detecting two low (~ 1 keV) tempera- 
tures close together. More counts at that temperature appear to be 
required to detect two temperatures. 

In order to examine whether there is significant non-radial 
variation in the bins, the fitted parameters (with l-o" errors) to each 
bin as a function of the radius of the bin are plotted in Fig. |l2| The 
radius plotted is the mean radius of the pixels of the bin. The plot 
brings out the features in Fig. ^| There is some significant variation 
in metallicity between radii of 7-50 arcsec. Examining those bins 
in the east and those to the west, we find that most of the variation 
in the abundance is to the west (the opposite side to the plume). The 
metallicity to the west is scattered between 0.5 and 1.7 Z Q at radii 
of 7-50 arcsec. The profile to the east is much smoother, except for 
an obvious high abundance point (Fig. |u](d), roughly 7 arcsec from 
the core). The east-west difference in metallicity and temperature is 
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Figure 10. (a) Adaptively binned image of the cluster, counts are average values per arcsecond pixel, (b) Upper-temperature results (for 2-T fits) or the 1-T 
fit results. The diagonal boxes inside each bin indicate the bounds of the l-o" error of the temperature. The white spots are excluded point-sources, (c) Lower- 
temperature fit results for the 2-T fits. White regions are where 1-T fits were used, (d) Abundance results for fits, (e) Absorption results of fits, (f) Reduced x 2 
of the fits. 
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Figure 12. Distribution of the fitted values for the bins in Fig. [nj as a func- 
tion of radius. The radius of a bin is the mean displacement from the centre 
of the cluster, and the error bars indicate the maximum and minimum radius 
of each pixels in the bin. The errors on the value are the l-o" fitted errors 
also shown in Fig. [Hi Two points are plotted in the temperature plot for bins 
fitted with a 2-T model. Bins at a radius of less than 4 arcsec are excluded, 
due to an incorrect fitted model (see text). 



highlighted by Fig. |j which shows east-side and west-side profiles 
of the cluster. 

Given that the cluster contains a strong radio source (Taylor et 
al. 2001) there may be significant non-thermal emission in the X- 
ray band. Any non-thermal emission will increase the continuum 
of the X-ray emission and lower the apparent metallicity. We tested 
for the presence of non-thermal emission by fitting the spectrum 
in each bin with a 1-T plus power-law model. We then used an 
F-test to decide whether adding the power-law component made a 
significant improvement to the fit (using the same criterion as § }4-.l| ) 
for each bin. 

We show in Fie. [l3]the areas where the power-law component 
is suggested by the F-test. Interestingly the morphology of the ra- 
dio source approximately matches those areas in the centre where 




n 1 r 

12h4Sra54s Kh+SmSls 12h48m48s 



Figure 13. Regions in the centre of the cluster where a power-law plus 1-T 
is an improvement over a 1-T fit. White areas indicate where the F-test gives 
no significant improvement over a 1-T fit. The other regions are shaded to 
show the photon index of the best-fitting power-law. The photon index was 
con trained to lie between 1.5 and 2.5. The contour marks where the radio 
emission is 0.2 mJy beam - 1 at 4.8 GHz, taken from the data of Taylor et al. 
(2001). 



the power-law component is a significant improvement. We also 
find that those areas 1 arcmin south of the core where the power-law 
is preferred match the low-intensity regions of the radio source (see 
Taylor et al. 2001 for detailed maps of the radio source morphol- 
ogy). However, only a small fraction of the emission from the core 
of the cluster appears to have a significant non-thermal component, 
and the peak in the abundance still exists if we exclude those re- 
gions where it may be important. Therefore the peak in abundance 
15 kpc from the core is not an effect of non-thermal emission. 



5 DISCUSSION 

The most interesting conclusion of our analysis is that the metallic- 
ity of the cluster peaks at ~ 1.3 Z Q (or 1.8 Z Q with a cooling flow 
model) at a radius of ~ 15 kpc, declining to 0.4 Z e at the centre. 
Previous observations (Fukazawa et al. 1994; Ikebe et al. 1998) of 
the cluster did not have high enough spatial resolution to resolve 
this feature, showing only an increasing gradient in metallicity to 
the centre. Bohringer et al. (2001) found using XMM-Newton data 
that the abundances in the Virgo Cluster peak ~ 1 arcmin from the 
centre. Molendi & Gastaldello (2001) later saw this to be an ef- 
fect of using a 1-T fit where more temperature components were 
required; within 1 arcmin from the centre the actual metal profile is 
flat. 

The abundance gradient in Centaurus was modelled by 
Reisenegger, Miralda-Escude & Waxman (1996), but no decrease 
was expected. We still get a peak in abundance just fitting spectra 
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Figure 14. Radial metallicity plot generated by fitting spectra between 3 to 
7 keV with a 1 -T model (above 3 keV only one temperature component is 
necessary). Each radial bin contains ~ 20000 counts from 0.5-7 keV. 



between 3 and 7 keV where Fe K emission occurs (Fig. |l4| ), so the 
peak is not an artefact of the Fe L complex. Resonant absorption is 
a possible mechanism for the metallicity decline in the core of the 
cluster. However, since we see the abundance peak in the Fe K lines 
too, resonant absorption is unlikely. 

In the future, it will be important to study the chemical com- 
position of the plasma, and handle projection effects fully. We also 
would like to understand why the statistical quality of the fits in to 
the spectra in the centre are poor, and how to improve them. These 
effects are probably due to the high metallicity of the core: a feature 
which is particular to this cluster. We would like to draw attention 
to the fact that there are differences in the qualities of fits using dif- 
ferent plasma codes. The fits to the adaptively-binned image spectra 
(Fig. UQKc)) suggest that the 2-T fits are required in specific regions, 
rather than at specific radii. 

Using a CF model produces a radial abundance profile with a 
similar shape, but with a higher peak metallicity, than that found 
using the 2-T model, indicating that the peak is not a model- 
dependent effect. The upper-temperature of the CF fit also matches 
the upper-temperature of the 2-T fit. In general, excluding the very 
central bins, the CF model gives a smaller reduced % 2 than the 2-T 
model. 

It is difficult to conclude in general which model gives the 
better fits to the regions of the cluster. When the data are binned 
in annuli, a standard cooling flow model (CF) is a significant im- 
provement (measured using an F-test with a maximum probabil- 
ity of 0.01) over a 1-T fit for nearly every annulus. 2-T and CFMT 
models are only a significant improvement over a 1-T model within 
about 60 arcsec of the cluster core. A 2-T model or CFMT model 
are a significant improvement to the CF model within the inner 
~ 7 arcsec. A CF model probably fits the spectra of the annuli bet- 
ter as the cluster is not radially symmetric, so there are a number of 
temperatures at each radius, which the CF model includes. 

The comparison gives different conclusions when fitting the 
regions in the adaptively binned image. A 2-T, CF or a CFMT 
model is a significant improvement (with the same F-test criterion) 
over the 1-T model inside much of the inner 40-50 arcsec. The area 
in the central arcmin for which the CF model is an improvement 
over the 1-T model is roughly 2/3 that for which the CFMT and 
2-T models are an improvement. The CFMT and 2-T models show 
a significant improvement in the fit over the CF model for areas 
about 50 arcsec away for the nucleus, and some of the area inside. 



Fig. |ll] shows that the radiative cooling time of the gas is less than 
a Gyr within those radii. Outside the core, there are some areas 
which show a significant improvement using a CF, 2-T or CFMT 
model over a 1-T model. It may be that there are local variations 
in temperature in these regions, or this could be a statistical effect. 
Using a larger bin size, with the corresponding increase in counts 
in the spectra, still results in areas outside the core which show a 
significant improvement using a 2-T or CFMT model over a 1-T 
model, although using a larger bin size includes more gas, which 
is probably at different temperatures. Note that if there is a cooling 
flow, then the models used here do not incorporate the gravitational 
work done on the gas in the flow. Therefore no firm statement can 
be made on the level of the mass deposition rate. The data support 
a rate of several tens M G yr~ 1 . 

It is difficult to speculate on the cause of the plume-like fea- 
ture in the core of the cluster. The galaxy velocity data cannot rule 
out a 100-200 km s" 1 line of sight motion of NGC 4696 in Cen 30 
similar to that of the cD in A1795 (Fabian et al. 2001b). The real 
velocity could of course be larger and there may be significant mo- 
tion of the intracluster gas within its potential well. The plume does, 
however, touch the east-most extent of the radio lobes (Taylor et al. 
2001), indicating there may be some interaction between the lobes 
and the plume. However, the direction of the plume is not aligned 
with the radio lobes. The plume has a higher metallicity than the 
core, suggesting the cooling flow is depositing enriched material in 
this area. 

Sparks, Macchetto & Golombek (1989) and de Jong et al. 
(1990) have proposed that the cooler central gas around NGC 4696 
originates from an infallen dust-rich dwarf galaxy, which also ac- 
counts for the dust lane. The fact that the metallicity of the plume 
matches its surroundings argues against this stripping hypothesis 
(as mentioned in §3.2). Moreover, cool gas is found in the majority 
of clusters. 

As in M87, the plume could be due to a radio bubble drag- 
ging out cooler gas (Bohringer et al. 1995; Churazov et al. 2001). 
The 'cold front' at the end of the plume could in fact be the edge 
of the bubble. However, we have shown that the metallicity of the 
plume is similar to the surrounding material, and higher than the 
central gas, which argues against that gas having originated close 
to the central radio source. Of course, the drop in metallicity within 
the inner 20 kpc (60 arcsec) could be due to bubbles transporting 
higher metallicity gas to larger radii (Churazov et al. 2001 ; Briiggen 
& Kaiser 2001). However, it is not then clear where the present 
lower-abundance gas has come from. More detailed models of the 
convection and exchange of gas are required, but at first inspection, 
such models are unlikely to reproduce the abundance profile seen 
here. 

Our main conclusion is that the steep abundance gradient in 
this cluster does not persist within 15 kpc, but decreases inward 
to a value similar to that seen at large radii. This may be a fea- 
ture of radiative cooling of the gas, if composed of small metal-rich 
clumps in a metal-poor surrounding environment (Morris & Fabian, 
in preparation). The presence of an abundance gradient means that 
there has not been significant widespread convection or mixing ex- 
tending beyond 15 kpc driven by the nucleus of NGC 4696. The 
spectra within ~ 60 kpc also indicate that the temperature structure 
of the gas has at least two phases, both decreasing towards the cen- 
tre. The present data are not capable of testing for the presence of 
further phases. 
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